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■ In the light of recent experimental data from the CLEO Collaboration we 

Q_i' study the decays of B mesons to a pair of pseudoscalar (P) mesons, and a 
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vector iy) meson and a pseudoscalar meson, in the framework of factorization. 
In order to obtain the best fit for the recent CLEO data, we critically examine 



^ ■ the values of several input parameters to which the predictions are sensitive. 



These input parameters are the form factors, the strange quark mass, ^ = 
\/N(. {N(. is the effective number of color), the CKM matrix elements and in 
particular, the weak phase 7. It is possible to give a satisfactory account of 
the recent experimental results in S — > PP and VP decays, with constrained 
values of a single ^. We identify the decay modes in which CP asymmetries 
are expected to be large. 
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I. INTRODUCTION 



The CLEO Collaboration has recently reported new experimental results on branch- 
ing ratios (BRs) of a number of exclusive decay modes where B decays into a pair of pseu- 
doscalars (P), a vector {V) and a pseudoscalar meson, or a pair of vector mesons. Several 
decay modes have been observed for the first time, such as B ^ vr+vr", K^ir^, ujtt^, p^vr^, 
p^TT^, K*^7t'^, and K*-)]. Improved new bounds have been put on the branching ratios for 
various modes, such as -B K^tt^, K^tt^, K^tt^, uK, uh^, and Kt]'. A search for CP 
asymmetries in S — K^tt^, K^tt^, K^n^, Kr]', and utt^ has also been performed. 

Recently, two works have been done to explain the recent CLEO results for B decays: one 
of them is based on the flavor SU(3) symmetry 0, and the other one involves the framework 
of factorization In the latter work, the existence of two different effective numbers of 
color, N^^{LL) and N^^{LR), is essential, and their favored values are N^^{LL) = 2 and 
N^^{LR) = 6 to explain the experimental data, except for B^ K*^7r^ and K^tt^. 

In the light of the recent CLEO data, in this work, we will analyze both B PP and 
B —>■ VP with a single parameter ^ in the framework of generalized factorization, in order 
to find satisfactory explanation compatible with all the recent experimental results. Our 
approach will be different from that of Ref. 0. Similar to the case oi B ^ D (i.e., heavy 

heavy) decays, we will assume only one universal ^ = 1/Nc (A^c is the effective number 
of color) in the analysis of both B PP and B VP {P and V are light mesons). It 
has been shown that in B ^ D decays a single ^ can satisfactorily explain experimental 
data for both B — > PP (such as Dn) and B —>■ VP (such as Dp) |^. In order to achieve 
this goal, the values of all the input parameters, e.g., the form factors, the strange quark 
mass, Nc, the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements etc., will be carefully 
examined and constraints on these parameters will be investigated. We shall see that it 
is indeed possible to account satisfactorily for all the recent experimental data within our 
framework. In addition, we will discuss CP asymmetries in the B decays and identify the 
decay modes where the CP asymmetries are possibly large. In our previous works P,p!0 



we 



have shown that the predictions for B PP and VP modes are sensitive to several input 
parameters, such as the form factors, the QCD scale, the parameter ^ = l/N^, the CKM 
matrix elements, and the light quark masses, in particular the strange quark mass. With the 
improved recent data, the results (e.g., ^ ~ 0) obtained in the previous works are unlikely to 
be compatible with the experimental results. Our main emphasis was to explain the large 
branching ratio for B — > Krj' within factorization. In order to explain the large branching 



ratio for B — > Krj', different assumptions have been proposed, e.g., large form factors |]TT 



the QCD anomaly effect |1T2| , [13| , high charm content in rj' ||I^-[T6|, a new mechanism in the 



Standard Model |T^, or new physics like supersymmetry without R-parity []TB[. In this work 
we will assume that a certain (unknown) mechanism is (at least in part) responsible for the 
large branching ratio for B Krj' and hence will not consider the decay modes having rf^'^ 
in the final states. 

We organize this work as follows. In Sec. II we discuss the effective Hamiltonian and 
obtain the effective Wilson coefficients at the scale for both 6 — >■ s and h ^ d transitions. 
In Sec. Ill we describe the parametrization of the matrix elements and define the decay 
constants and the form factors. In Sees. IV and V the two body decays B —>■ PP and 
B VP are analyzed in the framework of factorization. The CP asymmetries in the B 
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decay modes are also discussed. Finally, in Sec. VI our results are summarized. 



II. DETERMINATION OF THE EFFECTIVE WILSON COEFFICIENTS 

The effective weak Hamiltonian for hadronic B decays can be written as 



- 71" 



i=3 



+ H.C., (1) 

where Oj's are defined as 

0( = q^^Lfh^Lh, oi = qalf^Lffsfp^Lb^, 
03(5) = q'J^,LbT,q'y''L{R)q', ©4(6) = qal^lLhp'Lq'p'^^L{R)q'^, 

3 3 
^7(9) = -^Q7f,LbJ:eg>q'j''R{L)q', Os(w) = -ga7^L6^Seg/g^7''/2(L)g;^, 

On = J^m.qa^^RTabG^^'' , O^ = -^m.qa.^RbF^'' , (2) 

where L{R) = (1 =1= 75)/2, / can be m or c quark, q can be (i or s quark, and q' is summed 
over u, d, s, and c quarks, a and /5 are the color indices. is the SU(3) generator with 
the normalization Tr (T'T'') = 6°'^/2. G'^'^ and F^'^ are the gluon and photon field strength. 
q's are the Wilson coefficients (WC's). Oi and O2 are the tree level and QCD corrected 
operators. O^^q are the gluon induced strong penguin operators. Oy.io are the electroweak 
penguin operators due to 7 and Z exchange, and the box diagrams at loop level. In this 
work we shall take into account the chromomagnetic operator On, but neglect the extremely 
small contribution from O12. The dipole contribution is in general quite small, and is of the 
order of 10% for penguin dominated modes. For all the other modes it can be neglected. 
The initial values of the WC's are derived from the matching condition at the mw scale. 



However we need to renormalize them 19,20 when we use these coefficients at the mi, scale. 




We will use the effective values of WC's at the scale = rrib. It has been shown that in 
B — > PP case there is very little fj, dependence in the final states [^. 

We obtain the Cj(/i)'s by solving the following renormalization group equation: 

CK^^^/) = f^C{t,a,{fi),a,), (3) 

where t = ln(M^//i^) and C is the column vector that consists of (cj)'s. The beta and the 
gamma are given by 

where aem is the electromagnetic coupling and ^/ is the number of active quark flavors. 
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The anomalous-dimension matrices and 7^°^ determine the leading log corrections 
and they are renormalization scheme independent. The next to leading order corrections 
which are determined by 7^'' and 7^^-* are renormalization scheme dependent. The 7's have 
been determined in Refs. |T^PU| . 

We can express C(/i) (where fi lies between Mw and nib) in terms of the initial conditions 
for the evolution equations 



= U{fx,Mw)CiMw). 



(5) 



C{Mw)^s are obtained from matching the full theory to the effective theory at the Mw scale 
p0| , |2l[| . The WC's so far obtained are renormalization scheme dependent. In order to make 
them scheme independent we need to use a suitable matrix T 
/i = rrtfe are given by 



The WC's at the scale 



C{ix) = TUimb,Mw)C{Mw). 



The matrix T is given by 



.Op 



"Arc " 47r 



(6) 



(7) 



where f depends on the number of up-type quarks and the down type quarks, respectively. 
The r's are given in Ref. [^. In order to determine the coefficients at the scale /i < mi,, 
we need to use the matching of the evolutions between the scales larger and smaller than 
the threshold. In that case, in the expression for T, we need to use 6r instead of f, where 
= ru,d — ru,d-i (where u and d are the number of up type quarks and the number of 
down type quarks, respectively). The matrix elements [O'^s) are also needed to have one 
loop correction. The procedure is to write the one loop matrix element in terms of the tree 



level matrix element and to generate the effective Wilson coefficients p2| 



E 



Ci(/i) 



tree 



(8) 



where 



Ps 

P. 



-C2 



cq 



(3ci + C2] 



^eff 

C5 
eff 

Ce 

eff 
Cj 

eff 



10 



( 



Cl \ 

C3 - Ps/3 

C4 + Ps 
C5 - PsI'i 

ce + Ps 

Cj + Pe 

Cg + Pe 
ClO / 



+ G{mc,fi,q^ 
10 



+ 



VubV, 



uq 



VtbVtl 



cq 
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+ G{mc,fx,q ) + 



(9) 



10 

VubV. 



10 



(10) 



uq 



VtbVtl 
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+ G{mc,iJ,q 
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Here Vij are the elements of the CKM matrix, rric is the charm quark mass and is the 
up quark mass. The function G{m, yU, g^) is given by 



pi J 
G{m, jjL.q^) = A / dxx(l — a;)ln- 



where q is the gluon momenta in the penguin diagram ||l5| , ^ , p6| . In the numerical calcu- 
lation, we will use q^ = ml/2 which represents the average value and the full expressions 
for Ps e- In Table I we show the values of the effective Wilson coefficients at the scale rrib 
and mf,/2 for the process b — > sqq. Values for b — >■ dqq can be similarly obtained. These 
coefficients are scheme independent and gauge invariant. 

Since the color octet contribution is neglected in factorization approximation, we keep ^ = 
1/Nc as a variable. As an example, let us consider the u quark contribution of the operator 
O5 to the process B~ n^K'. There are two configurations: (part 1)=< K~\s'y^{l — 
75)6|i?~ >< 7r°|u7^(l + 75)u|0 >. For (part 2), we need to Fierz-transform the operator O5 
to -2{{1/N,)[uil - 75)&][s(l +75)m)] + (1/2)[m(1 - 75)Ae6] [s(l + 75)A,m]}, where is the 
color matrix. Only the (part 1) contributes in the factorization approximation. In order to 
account for this color octet term (part 2), one needs to make ^ a free parameter. 



III. HADRONIC MATRIX ELEMENTS IN FACTORIZATION 

APPROXIMATION 



The generalized factorization approximation has been quite successfully used in two body 
D decays as well as B ^ D decays [^]. The method includes color octet nonfactorizable 
contribution by treating ^ = l/A'c as an adjustable parameter |]7|j9| JT0|JT3|JT5|JT6|j23[| . In this 
work one of our goals is to establish the range of value of a single ^ for the best fit in both 
B — > PP and VP decays, where P and V are all light mesons such as tt, K^*\ p, u, and 0. 

Let us describe the parameterizations of the matrix elements and the form factors in the 
case of -B ^ PP and VP decays. 



{Pip')\V,\Bip)) 
{V{e,p')\{V,-A,)\B{p)) 



2 



Fiiq') + 



(12) 



TUB + my 

- [mB + my 
e*-q) 



p'^V{q^ 



AM' 



+ 



mB + my 



(P + P)/. 



A2{q') 



- (e ■ (l)^r%AQ{q 



and the decay constants, fp and /y, are given by 

m^\P{p)) = tfpp,, {0\V^\V{e,p)) = ifymye,, 



(13) 



where P, V , V^, and denote a pseudoscalar meson, a vector meson, a vector current, and 
an axial-vector current, respectively. p{p') and mp{my) are the momentum of B meson (P 
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or V) and the mass of P{V), respectively, q is given hj q = p — p' and is the polarization 
vector of V. Note that Fi(0) = Fo(0) and we set F{q^ = ml) = F{q^ = 0), since these 
form factors are assumed to be pole dominated by mesons at scale m^. Among all the form 
factors in the {V{e,p')\(y^ — A^)|-B(p)) matrix element, only Aq survives when we calculate 
the full B —>■ VP decay amplitude. The Aq is related to Ai and A2: 



Ao{0) 



niB + mv . , . niB 
-Ai[U) - 



my 



2m\ 



2m\ 



^2(0). 



(14) 



For our numerical calculations we use the following values of the decay constants (in 
MeV) i|igj2|: 



U = 132, fK = 162, /, = 215, = 215, = 225, = 237. 



(15) 



For the values of the form factors, two different sets, based on the Bauer-Stech-Wirbel (BSW) 
model |]2^ and light-cone QCD sum rule analysis |2^, have been often used in literature ||26|| . 
In our analysis, to achieve the best fit, we treat the form factors as parameters in a certain 
range of values which are reasonably consistent with the frequently used values. We first 
start with the following values of form factors : Fq^'^ = 0.33, Fq^^ = 0.38 obtained in the 
BSW model, and Aq^'' = Aq^'^ = 0.4. Then, to find the single ^ consistently explaining 
the experimental data, we vary the values of the form factors in a certain range as well as 
other parameters, such as the CKM weak phase 7 and the strange quark mass m^. 



IV. B DECAYS INTO TWO PSEUDOSCALARS 

Here we consider the decays B —>■ titt and Kn. The CLEO collaboration has made 
first observations of the decay modes B tt+tt^ and K^tt^, and has presented improved 
measurement of BRs for B —>■ K^n^ and B^ K^tt"^ and B"^ K'^ix^ as follows : 



B{B^ - 


-i> TT^TT ) 


= (4.3^i:^±0.5) X 10-^ 


BiB"" - 




= (17.2+1^ ± 1.2) X 10"^ 


B{B^ - 




= X 10^^ 


B{B^ - 




= (18.21^-^ ±1.6) X 10"^ 


B{B^ - 




= (14.6^^:?+^:^) X 10-^ 



(16) 



In Figs. 1—5, we plot the BRs averaged over particle- antiparticle decays for the modes 
7r+7r~, K^n^, K^ir^, K^ir^, and K^ir^ in final states as a function of ^ = l/N^ for /z = m^. 
In these figures, we use five different kinds of lines corresponding to different values of 
parameters, in addition to thick solid lines representing the experimental upper and lower 
bounds. The solid line corresponds to the case of choosing the form factors -Fq^^'^(O) = 0.33 
and F^^^{0) = 0.38 based on the BSW model, 7 = Arg(V;*J = 60°, m,(m;,) = 106 MeV, 
Vcb = 0.040, iKife/Vcbl = 0.087, and \Vtd\ = 0.004. The short dashed line then corresponds 
to the case of choosing 7 = 110" and the other parameters are the same as in the solid 
line case. Similarly, the dot-dash-dot line corresponds to the case of choosing ms{mb) = 85 
MeV, while the values of other parameters are the same as in the solid line case. For the 
dot-dashed and long dashed lines, we choose msirrib) = 106 MeV and 85 MeV, respectively, 
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with smaller values of form factors F(f ^^(0) = 0.26 and F(f^^(0) = 0.29, and 7 = 110°. 
Thus, by comparing each line with the solid line, one can easily see how the decay rate for 
any mode chang particular parameter, such as 7 or m^, changes. Note that in Fig. 1 

for — i> vr+vr^ decay, the dot-dash-dot and long dashed lines are identical to the solid and 
dot-dashed lines, respectively, since the amplitude for this decay mode does not depend on 
rris. Similarly, in Fig. 2 for 5^ K^tc^ decay, the short dashed line is identical to the solid 
line since the amplitude for this mode receives contribution from the penguin diagram only 
and is independent of 7. 

The decay rate for —>■ tt+tt^ is proportional to and is sensitive to the value 

of the form factor. In Fig. 1, one can see that the BR decreases, as the value of F(^^'^ 
decreases and/or the value of 7 increases. In order to fit the experimental upper limit on 
the BR for this mode, a smaller F^^'^ and a larger 7 (dot-dashed line, identical to the long 
dashed line) are favored. For F^'^'^{0) = 0.26 and 7 = 110°, the values of ^ ^ 0.44 are 
allowed. The rate for B^ K^n^ is also proportional to and is sensitive to 

as well. Figure 2 shows that the BR increases, as the value of F^^^'" increases and/or the 
value of rris decreases. In order to find a solution consistent with B^ tt+tt^, F^^'^ should 
not be too small and smaller rris is favored (long dashed line). For Fq^^'^(O) = 0.26 and 
^s{f^b) = 85 MeV (long dashed line), the allowed values of ^ are ^ ^ 0.45. However, for 
larger ms{mi,) = 106 MeV (dot-dashed line), the allowed values of ^ are smaller ^ < 0.15 
and these values of ^ are not allowed by the dot-dashed line for 5° tt+tt^. In Fig. 3, we 
plot the BR for 5° K'^tt^ averaged over particle-antiparticle decays as a function of ^ for 
fi = nib. This decay mode is sensitive to Fq^'^, 7 and m^. In the case of the long dashed 
line, the allowed values of ^ are ^ ^ 0.6, which are consistent with those in 5° tt+tt" and 
B~^ K^ir^. Similarly, in Figs. 4 and 5, we plot the BRs for K^ir^ and K^tt'^ in the final 
states as a function of ^ for /i = m^. These modes depend on Fq^'^, Fq^^, rrig, and 7. 
In the long dashed line case, ^ ^ 0.88 and ^ 0.47 are allowed for B K^tt^ and K^n^, 
respectively, and are consistent with those in the above decays (Figs. 1—3). Therefore, we 
conclude that the long dashed lines (the dot-dashed line in the case of i? ^ tt+tt") in Figs. 
1—5 represent the possible solution compatible with all the experimental limits on the BRs 
for decay modes tt'^7c~, K^n^, K^tt^, K^n^, and K^tc^ in the final states. The values of 
C, allowed by the data are those near C, ~ 0.45. (In fact, we shall see that the long dashed 
line can also consistently explain all the experimental data for B VP considered in next 
section.) Note that the modes B — >■ tt+tt", K'^ir^, and i^°7r° provide strong constraints on 
the values of ^ to satisfy the data. In particular, the data on the modes B — >■ 7r"'"'7r~ and 
K^n^ put very tight limits on the allowed values of ^. Hence, an improved measurement of 
these modes will be important in testing the framework of factorization. 

The experimental bounds on the BRs for decays B —>■ PP provide constraints on the 
parameters. The favored values of the parameters are 

Fo^^"(0) = 0.26, F(f^^(0) = 0.29, (17) 

7 ^ 110°, m,(mb) = 85 MeV, 

Kb = 0.040, iKb/Kfel = 0.087, 11^^,1= 0.004. 

Current best estimates for CKM matrix elements are Vet = 0.0381 ± 0.0021 and = 
0.085 ± 0.019 [^. The CLEO has recently made the first determination of the value of 
7 = 113°^230 by any method other than the unitarity triangle construction []5|, p8| , p9| . The 
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favored values for the CKM matrix elements in our analysis above are chosen to get the best 
fit for the experimental limits on the BRs for the decay processes B — > PP and VP. We 
find that, if \Vub\ increases, the rates for tt+tt" and uh {h = vr, i^) increase, while the rate 
for K'^n^ decreases. Also if \Vtd\ increases, then the rates for tt^tt^ and K'^tt^ increase and 
the rate for ujh decreases. 

In Table II, we present the BRs and the CP asymmetries for B PP decays at a 
representative value of ^ = 0.45. (We shall see that the values of ^ near ^ ^ 0.45 are favored 
to fit all the data.) Available experimental values are also presented. The BRs for all the 
modes are compatible with the present experimental data. The CP asymmetry, Acp, is 
defined by 

_ B(i, /) - B{b ^ /) 

where h and / denote h quark and a generic final state, respectively. The recent CLEO search 
for CP asymmetries in i? ^ Ktt decays has found : —0.70 < Acp{B^ —>■ K^n^) < 0.16, 
-0.35 < Acp{B^ K^TT^) < 0.27 at 90% confidence level (C.L.). The expected CP 
asymmetries in B ^ PP decays are generally small and range from —11% to 0. 

The BRs in this analysis have been evaluated at the scale mi,. However, if we had chosen 
the QCD scale /i = mb/2, the result would not change much [Q. We also see that the favored 
values of parameters involve a lighter strange quark mass. This, however, is in accordance 
with the latest trend of lattice results |3^. The ratios of the quark masses are much better 



known than the individual masses. For example [pi|], 

TTl Tfl 

— = 0.553 ±0.043; — = 18.9 ± 0.8 (19) 

rud rud 

The strange quark mass is in considerable doubt: i.e., QCD sum rules give ms{l GeV) = 
(175 ± 25) MeV and lattice gauge theory gives m^(2 GeV) = (100 ± 20 ± 10) MeV in the 
quenched lattice calculation In this analysis we have varied from 150 to 116 MeV 
at 1 GeV scale. We see that the of 116 MeV gives rise to the best fit. We have used the 
quark masses at the mi, scale. The magnitude of reduces from 150 to 106 MeV and from 
116 to 85 MeV at the mb scale through 3 loop QCD and 1 loop QED RGEs. The magnitudes 
of the other quark masses (m^ and m^) also depend on the strange quark mass. Satisfying 
the constraints from Eq. (p!9D, the values of md we have used are 5.9 MeV (corresponding 
to ms=150 MeV) and 4.7 MeV (corresponding to ms=116 MeV) at the scale. Similarly, 
the values of m„ we have used are 3.4 MeV (corresponding to ms=150 MeV) and 2.8 MeV 
(corresponding to m<j=116 MeV) at the scale. 

The difference in affects the BRs for the B — > t^K modes, since the amplitudes 
contain a factor such as X = M]^/{{mb + mu){ms + m„)). On the other hand, changes in 
md or m„ do not affect the BRs significantly. In the case of \^S\ = 1 decays, m„ or can 
always be neglected compared to or in the factor X. In the case of = decays, 
the tree contribution is large compared to the penguin contribution and the factor similar 
to X appears only in the penguin term. Hence the effect is small. For example, the BR for 
B^ TT+TT^ changes from 6.00 x 10~^ to 6.10 x 10~^ at = 0.45, when one changes the 
from 5.9 MeV to 4.7 MeV (changing the m^ mass as well). Similar results hold true also in 
the B — > VP case. 
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V. B DECAYS INTO A VECTOR AND A PSEUDOSCALAR 



We now analyze the decay processes B —>■ VP which include B —>■ uJ7r{K), p7t{K), 
(f)'7i{K), and K*n{K). The recent measurement at CLEO has yielded the following bounds 
i : 

(20) 



B{B^ 




= (11-311^ 


± 1.4) 


X 


10 


B{B^ 




= (14.311^ 


±2.0) 


X 


10 


B{B^ 




= i^OA+li 


± 2.1) 


X 


10 


BiB" - 






±4.2) 


X 


10 


BiB" ^ 




= mtitt) 


X lo~^ 





and 

B{B^ ujK^) < 7.9 X 10~^ (21) 
B{B^ ^ 0K±) = (GAtlltli) X 10-^ 

where denotes tt^ or K^, and the BR for B^ p^vr^ is the sum of the BRs for 
B° p+TT" and 5° p~7r+. Note that the above BR for 5° K^^n"^ still involves a 
large error. 

As in the case of -B ^ PP decays, in Figs. 6—11, we plot the BRs averaged over 
particle-antiparticle decays for the modes B —>■ utt"^, uh'^, uK"^, p^ii^ ., (pK^, and i^*^7r'^ 
as a function of ^ for fi = mf,. Six different kinds of lines are used, corresponding to different 
values of parameters. The definitions of the (five) lines are the same as those in 5 — PP 
case, except that the form factors Aq{0) = Aq^'^{0) = Aq^'^{0) = 0.4 are now added. The 
dotted line is newly introduced, which corresponds to Aq{0) = 0.36 with the same values of 
other parameters as those in the solid line case. Thus, a comparison of the dotted line with 
the solid line shows how the BR for a particular mode changes as Aq changes. 

In Fig. 6, we present the plot of the BR for B^ — > cuvr^ as a function of ^. This decay 
mode receives the dominant contribution from the tree diagram and is sensitive to the form 
factors Aq^'^, F^^'^, and the weak phase 7. The long dashed line and the dot-dash-dot 
line are identical to the dot-dashed line and the solid line, respectively, since the rate for 
this mode does not depend on m^. All the lines are well within the experimental limits for 
values of ^ in a broad region. In the case of the long dashed line, the allowed values of ^ are 
0.3 ^ ^ ^ 0.75. The recent CLEO search for CP asymmetry in B^ uir^ decay has found 
: -0.80 < AcpiB^ ujTT^) < 0.12 at 90% C.L. We find that the expected CP asymmetry 
in this mode is —11% for a representative value of ^ = 0.45 (we shall see below that the 
values of ^ near ^ ^ 0.45 are the favored values for the best fit). 

The plot of the BR for B^ ujK^ as a function of ^ is shown in Fig. 7. The rate for 
this process depends on Aq^'^, F^^^ , 7, and m^. The previous experimental result from 
CLEO for this decay mode ^ showed the large BR of BiB^ uK^) = (15l^ ±2) x 10" 



but in the recent CLEO report [0 the statistical significance for this mode is only 2. la and 
the upper limit for the BR at 90% C.L. has been set as in Eq. (|2l|), which is much lower 



than the previous one. Thus, as can be seen in Fig. 7, the values of ^ in a broad region 
are compatible with the experimental upper limit. However, in our earlier work [|l^], only 
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smaller values of ^ ^ 0.05 or larger values of ^ > 0.6 were allowed to fit the previous data. 
The allowed values of ^ for the long dashed line are ^ ^ 0.67. At a representative value of 
^ = 0.45, the expected BR for this mode is B{B^ ojK^) = 1.33 x 10"^. 

Figure 8 shows the plot of the BR for uh^ as a function of ^, where /i is vr or 

K. To obtain the best fit for this mode, smaller values of Ff~*^, F^~*^ , and 7, and larger 
values of Aq^'^ and are favored. The previous CLEO measurement of the BR for this 
mode 15211 was B{B^ uoh^) = (25^7 ±3) x 10"^ but the recent value § has been reduced 
to B{B^ uh^) = (14.3l|^ ± 2.0) x 10"^. Thus, the values of ^ < 0.1 and 0.4 < ^ < 0.64 
for the long dashed line are compatible with the recent data in this mode, while the previous 
allowed values of ^ are ^ ~ and C, ^ 0.5. The expected BR for this mode at a representative 
value of ^ = 0.45 is B{B^ uh^) = 11.34 x 10"^ 

The case of B^ (pK^ is shown in Fig. 9. The previous CLEO result for this mode 
H was B{B^ (pK^) < 5 X 10"^ at 90 % C.L. but recently CLEO has announced the 



new data : B{B^ (pK^) = (6.4l|?l^;^) x 10"^, 0A^O) < 12 x 10-^ and the 

combined branching ratio B{B 4>K) = (6.2lf;gl°;7) x 10^^. (The Belle Collaboration has 
recently reported the branching ratio : B{B^ 4>K^) = (17.2j;5;J ± 1.8) x 10"®, which 
is very large and inconsistent with the CLEO data. To be consistent, in this analysis we use 
the recent CLEO data only. Future improved data for this mode are called for.) This decay 
is a pure penguin process and is sensitive to Ff^^, but independent of Aq, 7, and m^. A 
smaller F^^^ is favored for a better fit in this decay. But the decays B Kir disfavor 
too small values of F^{^'' and F^{^^ . We find that Fo^^^(O) ^ 0.26, F^^^{0) ^ 0.29 are 
favored. For the dot-dashed line, identical to the long dashed line in this mode, the values 
of 0.3 ^ ^ ^ 0.55 are compatible with the data. 

In Fig. 10, we present the cases of B ^ pn decays in which the tree contribution is 
dominant. In (a), the sum of the BRs for B^ p+7r" and B^ p"7r+ as a function of ^ is 
shown in order to compare with the recent CLEO data in Eq. (|20|) . The decay B^ p^vr" 
is sensitive to Ff^'^ and 7, while 5° p'Tr^ is sensitive to Aq^^ and 7. The lines are well 
within the experimental upper and lower limits for most values of C,- For the dot-dashed 
line, identical to the long dashed line, the allowed values of ^ are ^ > 0.08. The expected 
BR for this mode is B{B^ p^vr^) = 29.41 x 10~® at a representative value of = 0.45. 
The case of B^ p^ir^ is shown in (b). This process is sensitive to Aq^'', F^^'^, and 
7. For the dot-dashed line, identical to the long dashed line, the values of ^ ^ 0.34 are 
compatible with the experimental data. This can be compared with the case of B^ un^ 
shown in Fig. 6. These two modes are both tree-dominated and have similar values of 
masses, decay constants, and form factors. So it is expected that their BRs are not very 
different : the recent CLEO data for p^vr^ and lutt^ decays are (10.4l|l ± 2.1) x 10"® and 
(11.3tO ± 1-4) X 10"®, respectively, as in Eq. (^). The favored values in our analysis are 
^ p°7r±) = 10.06 X 10"® and un^) = 10.01 x 10"® at a representative value 

of ^ = 0.45, which are consistent with the recent data. 

Figure 11 shows the plot of B{B^ K*^^^) as a function of ^. This mode receives the 
dominant contribution from the penguin diagram. Our theoretical expectation of the BR 
for this decay is less than the experimental limits at la level (thick lines), but is greater 
than the lower limit at 2a level (gray line). Thus, within 2a range, our result is compatible 
with the data for this mode. Since the measurement of this decay still involves large error, 
an improvement in the experiment will be crucial to test the framework of our work. 
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In Tables III and IV, we present the BRs and the CP asymmetries for the decays B VP 
{AS = and \ AS\ = 1) at a representative value of ^ = 0.45. Available experimental results 
are also presented. All the theoretical values are compatible with the present experimental 
bounds. In particular, in some decay modes, the CP asymmetries are expected to be large. 
Among AS = decays, the CP asymmetry is expected to be relatively large in a few decay 
modes: (i) in p'^^r^, the expected CP asymmetry is —27% with the expected BR of 

0.65 X 10~^, (ii) in B^ ujtt^, the CP asymmetry is expected to be —18% with the expected 
BR of 0.032 X 10^^, (iii) in B^ ujtt^, the CP asymmetry is expected to be —11% with 
the expected BR of 10.01 x 10~^. Among \AS\ = 1 decays, there are several interesting 
modes: (i) the expected CP asymmetry in B^ ujK^ is —29% with the expected BR of 
0.15 X 10~^, (ii) the expected CP asymmetry in B^ uoK^ is —19% with the expected BR 
of 1.33 X 10~^, (iii) the CP asymmetry in 5° K*^tx^ is expected to be —14% with the 
expected BR of 6.27 x 10~^, (iv) the CP asymmetry in 5^ K*^tt^ is expected to be 15% 
with the expected BR of 3.56 x 10~^. 



VI. CONCLUSION 

Motivated by the recent CLEO data, we have analyzed charmless hadronic two body 
decays of B mesons B PP and B VP. In the framework of generalized factorization, 
we have carefully examined the values of several input parameters to which the predictions 
are sensitive. Those input parameters are the form factors, the strange quark mass, ^ = 
1/A^C5 the CKM matrix elements, and in particular, the weak phase 7. 

We have found that the experimental bounds on the BRs for the decay modes B tt+tt^, 
K^TT^, and K^n^ among B —>■ PP modes put strong constraints on the parameters. The 
constraints on parameters from the decays B^ — > uh^ among B —>■ VP are also strong and 
lead to the following favored values of the parameters for the best fit (in Figs. 1—11, the 
long dashed line (or the dot-dashed line when the long dashed line is absent) represents the 
case corresponding to the best fit): 

^ ^ 0.45 (22) 

F(f-"(0) = 0.26, F(f^^(0) = 0.29, 

A^^%0) = 0.4, A^^'^(O) = 0.4, 

7 ^ 110°, msirrib) = 85 MeV, 

Vcb = 0.040, \Vub/Vb\ = 0.087, \Vtd\ = 0.004. 

[In fact, a little smaller values of Ao{0) are also allowed.] It has been known that there 
exists the discrepancy in values of 7 extracted from the CKM-fitting at p — 77 plan [^] and 



from the analysis of hadronic decays of B mesons p9| . The value of 7 obtained from the 



each case is 7 = 60° ~ 80° from the CKM-fitting a.t p — rj plane, or 7 = 90° ~ 140° from 
the hadronic B decay analysis. In our analysis we find that 7 ~ 110° is favored to fit the 
given data. We have shown that the recent CLEO data in B —>■ PP and VP modes can be 
satisfactorily explained with ^ ^ 0.45, except for the BR of the decay mode -B° — > K*^7r^ 
at la level [at 2a level, our prediction for B{B^ —>■ K*^'n-^) is compatible with the data]. An 
improved measurement of the BR for this process will be crucial in testing the framework 
of factorization. We have also identified the decay modes where the CP asymmetries are 
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expected to be large, such as -B — > p°7r°, a;7r°, a;7r^ in AS" = decays, and B — > cuK'^, cuK^, 
K*^n°, p°K° in \AS\ = 1 decays. 
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TABLES 



TABLE L Effective Wilson coefficients for the b ^ s transition at the scales = and mb/2. 



WC's 




H = mb 


H = mb/2 


„eff 
Cl 

„efr 

C2 

„efr 

C3 

„efr 

C4 

„efr 

„efr 

^6 

„efr 

C7 

„efr 

^8 

„efr 

„efr 




1.149 

— u.ozuy 
0.02175 - 0.00414i 
-0.04906 - 0.01242i 
0.01560 + 0.00414i 
-0.06063 - 0.01242Z 
-0.000859 + 0.000073? 

0.001433 
-0.011487 + 0.000073i 
0.003174 


1.135 

— U.ZoZ 

0.0228718 + 0.004689i 
-0.051144 - 0.004689i 
-0.051144 - 0.004689? 
-0.0653549 - 0.0140673? 
0.00122773 + 0.00005724? 
-0.0000953211 
-0.0120155 + 0.0000572433? 
0.00218628 


TABLE IL The branching ratios (B) and the CP asymmetries (Acp) for B decay modes into 
two pseudoscalar mesons at a representative value of ^ = 0.45. Charged conjugate modes are 
implied. The experimental Acp's represent 90% confidence level intervals. 


Decay modes 


B (10-^ 


') Experimental B {10~^) 


Acp Experimental Acp 


B+ ^ ;r+7rO 

^ ^+^- 


4.38 
6.00 
0.11 


< 12.7 
4.3iJ:^ ±0.5 

< 5.7 




-0.11 



B+ ^ ii:07r+ 
50 ^ ^0^0 


10.14 
13.34 
15.90 
8.70 


11 f;+3.0+1.4 

18.2tto ± 1-6 
17.2111 ± 1.2 

14 «+5.9+2.4 


-0.061 [-0.70,0.16] 


-0.066 [-0.35,0.27] 
-0.01 
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TABLE III. The branching ratios {B) and the CP asymmetries (Acp) for B decay modes 
(AS = 0) into a vector and a pseudoscalar meson at a representative value of ^ = 0.45. Charged 
conjugate modes are imphed. The experimental branching ratio with ^ below is the sum of 
B° P'^tt" and p'~Tr'^. The experimental Acp represents a 90% confidence level interval. 



Decay modes 


B (10-6) 


Experimental B (10-^) 


Acp 


Experimental Acp 


B+ ujTr+ 


10.01 


11.3l|^ ± 1.4 


-0.11 


[-0.80,0.12] 


B+ ^ pOTT+ 


10.06 


10.4l|^ ± 2.1 


0.047 




B+ ^ p+ttO 


10.51 


< 43 


-0.039 




B+ K*^K+ 


0.043 


< 5.3 







B+ K*+K^ 


0.013 









B+ (f)Tr+ 


0.0028 


< 4.0 








5° ^ CUTT^ 


0.032 


< 5.5 


-0.18 


B^ ^ p+TT- 


15.39 


27.61? j ± 4.2t 


-0.072 


SO ^ p-7r+ 


14.02 




-0.015 


B^ ^ pO^o 


0.65 


< 5.5 


-0.27 




0.13 







B^ ^ <^7rO 


0.0014 


< 5.4 






TABLE IV. The branching ratios {B) and the CP asymmetries (Acp) for B decay modes 
(lA^I = 1) into a vector and a pseudoscalar meson at a representative value of ^ = 0.45. Charged 
conjugate modes are implied. The experimental branching ratio with * below is from Ref. [34]. 



Decay modes 


B (10-6) 


Experimental B (10 6) 


Acp 


B+ ojK+ 
B+ ^ pOK+ 

B+ p-K^ 

B+ K*^TT+ 
B+ K*+TT^ 

B+ (t>K+ 


1.33 
0.67 

1.20 
4.06 
3.56 
6.56 


< 7.9 

< 17 

< 48* 

< 16 

< 31 

f- 4+2.5+0.5 

O-"i-2.1-2.0 


-0.19 
0.13 




0.15 




5° p-K+ 
fiO ^ pOi^o 

B^ K*+TT- 
5° ^ cl>K^ 


0.15 
1.28 
0.082 
0.76 
6.27 
6.56 


< 21 

< 32 

< 27 

< 3.6 

99+8+4 
^^-6-5 

< 12 


-0.29 
0.10 
0.17 
-0.13 
-0.14 
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FIG. 1. Branching ratio for 



The solid line : Fo^^'^(0)=0.33, F^' 



The short dashed line 
The dot-dashed line : 



^0 



(0) 



-(0) 
'^(0)=0.26 



0.33, Fr^ 



^0 



as a function of .^(= -^). 
0.38, 7 = 60° and ms{mb)=106 MeV. 
^(0)=0.38, 7 = 110° and m^(mb)=106 MeV. 
^(0)=0.29, 7 = 110° and m,{mb)=106 MeV. 
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FIG. 2. Branching ratio for 



K^TT^ as a function of 

_ , =85 MeV. 

The dot-dash-dot line : F(f ^''(0)=0.33, Fo^^-^(0)=0.38, 7 = 60° and ms{mb)=85 MeV. 
The definitions for other hnes are the same as those in Fig. 1. 



The long dashed line : Fo^~'''(0)=0.26, F(f ^■^(0)=0.29, 7 = 110° and ms{mb) 
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FIG. 3. Branching ratio for K^ir^ as a function of -jy-). 

The definitions for the lines are the same as those in Figs. 1 and 2. 
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FIG. 4. Branching ratio for K^tt^ as a function of 

The definitions for the hnes are the same as those in Figs. 1 and 2. 
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FIG. 5. Branching ratio for — > K^tt'^ as a function of -^) 
The definitions for the hnes are the same as those in Figs. 1 and 2. 
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^ - 

0.2 0.4 0.6 0.8 1 

FIG. 6. Branching ratio for ujtt^ as a function of ^(= ■^). 

The sohd Hne : Fo^^'^(0)=0.33, F(f-"^(0)=0.38, 7 = 60°, Ao{0) = 0.4 and m^(mb)=106 MeV. 
The short dashed line : F(f ^'^(0)=0.33, F(f ^■^(0)=0.38, 7 = 110°, Ao{0) = 0.4 and m^(mb)=106 
MeV. 

The dotted line : F(f ^'^(0)=0.33, F^f ^■^(0)=0.38, 7 = 60°, ^o(O) = 0.36 and ms(mb)=106 MeV. 
The dot-dashed line : Fo^^'^(0)=0.26, Fo^^-^(0)=0.29, 7 = 110°, Ao{0) = 0.4 and m^(m,,)=106 
MeV. 
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b) 

FIG. 7. Branching ratio for — > loK^ as a function of .^(= ^j^). 
The long dashed hne : Fo^^'^(0)=0.26, Fo^^^(0)=0.29, 7 = 110°^ Ao(0) = 0.4 and m^(mb)=85 
MeV. 

The dot-dashed-dot hne : F(f ^'^(0)=0.33, Fo^^-f^(0)=0.38, 7 = 60°, ^o(O) = 0.4 and m,(mb)=85 
MeV. 

The definitions for other hnes are the same as those in Fig. 6. 
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FIG. 8. Branching ratio for — > w/i^ as a function of .^(= jf\ 
The definitions for the hnes are the same as those in Figs. 6 and 7. 
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FIG. 10. Branching ratios for i?^ p^vr^ and p^ir^ as a function of 

The definitions for the Unes are the same as those in Figs. 6 and 7. 



in 



o 

A 
ol 
CO 



2 . 5 



1 . 5 




0.2 



. 4 



. 6 



. 8 



FIG. 11. Branching ratio for K*^7r^ as a function of -^). 

The definitions for the Unes are the same as those in Figs. 6 and 7. 
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